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<^ • Abstract. Interstellar grains span a wide range of sizes from a few 

\ angstroms to a few micrometers. The presence of nanometer-sized or 

smaller particles in the interstellar medium is indicated directly by the 
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, interstellar far ultraviolet (UV) extinction, the ubiquitous 3.3, 6.2, 7.7, 
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8.6, and 11.3 //m polycyclic aromatic hydrocarbon (PAH) emission fea- 



^ , tures, the near and mid infrared broadband emission seen in the IRAS 12 

t::^- ; and 25 ^um bands and the COBE-DIRBE 3.5, 4.9, 12 and 25 ^m. bands, 

the 10-100 GHz Galactic foreground microwave emission, and indirectly 
by the heating of interstellar gas. For nanoparticles under interstellar 
conditions, UV/visible photon absorption is the dominant excitation pro- 
\^ I cess. With a heat capacity smaller than or comparable to the energy 

. of an energetic stellar photon, nanoparticles are subject to single-photon 

heating, followed by vibrational relaxation, photoionization, and pho- 
todestruction. With excited electrons spatially confined, semiconductor 
nanoparticles are expected to luminesce efficiently. This review focuses on 
the photophysics of nanoparticles with emphasis on the stochastic heat- 
ing and the vibrational excitation of interstellar PAH molecules, and the 
excitation of photoluminescence with special attention given to silicon 



O \ nanoparticles. 

to 
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1. Introduction: Historical Perspectives of Studies of Interstellar 
Grain Compositions and Sizes 

x: 

^ , When the existence of solid particles in interstellar space was first convincingly 

demonstrated by Trumpler (1930) through the discovery of color excess be- 
tween the photographic (with an effective wavelength Agg ~4300A) and visual 
(Aofr«^5500A) magnitudes, the grains were estimated to be ;^10 g, corre- 
sponding to radii of ^^20 A (Trumpler 1930). By the end of the 1930s, a 
extinction law (i.e., the interstellar extinction varied approximately inversely 
with wavelength A) in the wavelength range 1-3 ^m~^ had been well established 
(Hall 1937; Stebbins, Buffer, & Whitford 1939). At that time, metalhc grains 
with a dominant size of ~ 0.05 /im (Schalen 1936) or a power-law size distribu- 
tion dn{a) / da ^ a~^'^ in the size range 80A^a;^lcm (Greenstein 1938) were 
proposed to explain the A~^ extinction law, partly because meteoritic particles 
(predominantly metallic) and interstellar grains were then thought to have the 
same origin. 

In view of the spatial correlation between gas concentration and dust extinc- 
tion, Lindblad (1935) suggested that interstellar grains were formed by conden- 
sation from the interstellar gas through random accretion of gas atoms. However, 
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it was found later that in typical interstellar conditions, the Lindblad conden- 
sation theory would result in a complete disappearance of all condensable gases 
and the grains would grow to sizes (~ 10 fim) well beyond those which would 
account for the interstellar extinction (see Oort & van de Hulst 1946 and refer- 
ences therein). By introducing a grain destruction process caused by grain-grain 
collisions as a consequence of interstellar cloud encounters, Oort & van de Hulst 
(1946) further developed the interstellar condensation idea and led to the "dirty 
ice" model consisting of saturated molecules such as H2O. CH,i, and NH3 with 
an equilibrium size distribution which could be roughly approximated by a func- 
tional form dn (a) /da ~exp — 5 (a/0.5 /Ltm)^ and an average size of ~ 0.15 /xm. 

The discovery of interstellar polarization (Hall 1949; Hiltner 1949) cast 
doubts on the "dirty ice" model since ice grains are very inefficient polarizers. 
This stimulated Cayrel &: Schatzman (1954) to consider graphite grains as an 
interstellar dust component since aligned graphite particles would be a very 
efficient interstellar polarizer because of their strong optical anisotropy. The 
"dirty ice" model was also challenged by nondetection of the 3.1 /nm feature of 
H2O ice outside of molecular clouds (Danielson, Woolf, & Gaustad 1965; Knacke, 
Cudaback, & Gaustad 1969). This gave the late J. Mayo Greenberg (1922- 
2001) the incentive to perform the early experiments on the ultraviolet (UV) 
photoprocessing of low temperature mixtures of volatile molecules simulating 
the "original" dirty ice grains (Greenberg ct al. 1972) to understand how and 
why the predicted H2O was not clearly present. From such experiments was 
predicted a new component of interstellar dust in the form of complex organics 
known as "organic refractories" , containing a mixture of aliphatic and aromatic 
carbonaceous molecules (Greenberg et al. 2000). 

As an alternative to the interstellar condensation process, Hoyle &; Wickra- 
masinghe (1962) proposed that graphite particles with radii a few times 0.01 /xm 
could form in the atmospheres of cool N-type carbon stars, and subsequently 
be ejected into interstellar space by the stellar radiation pressure, and be re- 
sponsible for (part of) the interstellar extinction and polarization. The graphite 
proposal further gained its strength from the prominent 2175 A interstellar ex- 
tinction hump detected by the Aerobee rocket observation (Stecher 1965), which 
was generally attributed to small graphitic grains (Stecher & Donn 1965), al- 
though its exact nature still remains unidentified (Draine 1989, 2003a). 

Similarly, Kamijo (1963) first proposed that Si02 whose size is of the or- 
der of ~ 20 A, condensed in the atmospheres of cool M-type stars and expelled 
into interstellar space, could serve as condensation nuclei for the formation of 
"dirty ices". It was later shown by Gilman (1969) that grains around oxygen- 
rich cool giants are mainly silicates such as Al2Si03 and Mg2Si04. Interstellar 
silicates were first detected in emission in the Trapezium region of the Orion 
Nebula (Stein & Gillett 1969); in absorption toward the Galactic Center (Hack- 
well, Gehrz, Sz Woolf 1970), and toward the Becklin-Neugebauer object and the 
Kleinmann-Low Nebula (Gillett & Forrest 1973). Silicates are now known to be 
an ubiquitous interstellar dust component (see §1 in Li & Draine 2001a for a 
review). 

In the 1960s and early 1970s, the extension of the extinction curve toward 
the middle and far UV (A^3 fim^^) was made possible by rocket and satellite ob- 
servations, including the rocket-based photoelectric photometry at A = 2600 A 
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and 2200 A (Boggcss & Borgman 1964); the Aerobee Rocket spectrophotome- 
try at 1200 A ^ A ^3000 A (Stcchcr 1965); the Orbiting Astronomical Satellite 
(OAO-2) spectrophotometry at 1100 A ^ A ^ 3600 A (Bless & Savage 1972); and 
the Copernicus SatelHte spectrophotometry at 1000 A ^ A ^ 1200 A (York et al. 
1973). By 1973, the interstellar extinction curve had been determined over the 
whole wavelength range from 0.2 /Ltm~^ to 10 /xm~^. The fact that the extinction 
continues to increase in the far UV (e.g., see York et al. 1973) implies that no 
single grain type with either a single size or a continuous size distribution could 
account for the observed optical to far-UV interstellar extinction (Greenberg 
1973). This led to the abandonment of any one-component grain models and 
stimulated the emergence of various kinds of models consisting of multiple dust 
constituents, including silicate, silicon carbide, iron, iron oxide, graphite, dirty 
ice, solid H2, etc. 

In 1956 Piatt first suggested that very small grains or large molecules of 
less than 10 A in radius grown by random accretion from the interstellar gas 
could be responsible for the observed interstellar extinction and polarization. 
Piatt (1956) postulated these "Piatt" particles as quantum- mechanical particles 
containing many ions and free radicals with unfilled electronic energy bands. 
Donn (1968) further proposed that polycyclic aromatic hydrocarbon-like "Piatt 
particles" may be responsible for the UV interstellar extinction. 

Greenberg (1968) first pointed out that very small grains with a heat content 
smaller than or comparable to the energy of a single stellar photon, cannot be 
characterized by a steady-state temperature but rather are subject to substan- 
tial temporal fluctuations in temperature. Under interstellar conditions, grains 
larger than a few hundred angstroms in radius would attain an equilibrium tem- 
perature in the range of 15-25 K (Greenberg 1968, 1971; Mathis, Mezger, & 
Panagia 1983; Draine &: Lee 1984; Li & Draine 2001b). This temperature was 
long thought to be ~3K before van de Hulst (1949) made the first realistic 
estimation of 10-20 K for both metallic and dielectric grains. 

In a study of the scattering properties of interstellar dust (albedo and phase 
function) determined from the OAO-2 observations at 1500 A ^ A ^ 4250 A of 
the diffuse Galactic light (Witt & Lillie 1973), Witt (1973) first explicitly sug- 
gested a bi- modal size distribution for interstellar grains: large grains with radii 
;^2500 A would provide extinction in the visible region including scattering which 

is strongly forward directed, and small particles with radii ^250 A would domi- 
nate the UV region and contribute nearly isotopic scattering. 

The presence of "Piatt" particles in a dust cloud near M17 was first sug- 
gested by Andriesse (1978) based on an analysis of its IR spectral energy distri- 
bution which was shown to be characterized by a combination of widely different 
color temperatures (e.g., the far-IR and near-IR color temperatures are '~38K 
and ~150K, respectively) and that the 8-20 /xm emission spectrum is similar 
over a distance of ^ 2' through the source. Andriesse (1978) argued that the color 
temperature differences cannot be easily ascribed to a spatial variation of the 
dust temperature in this cloud, because the cloud profile appears to be indepen- 
dent of the wavelength; instead, "Piatt" particles with sizes of ~ 10 A exhibiting 
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temperature fluctuations up to ~ 150 K plus large grains with an equilibrium 
temperature of ~ 36K can explain the observed IR spectrum of M17. 

In a near-IR (1.25-4.8 /im) photometric and spectrophotometric study of 
three visual reflection nebulae NGC 7023, 2023, and 2068, Sellgren, Werner, & 
Dinerstein (1983) discovered that each nebula has extended near-IR emission 
consisting of emission features at 3.3 and 3.4 yUm and a smooth continuum which 
can be described by a greybody with a color temperature ^ 1000 K. They found 
that the emission spectrum (i.e., the color temperature for the continuum and 
the 3.3 /im feature) shows very little variation from source to source and within 
a given source with distance from the central star. Sellgren (1984) argued that 
this emission could not be explained by thermal emission from dust in radiative 
equilibrium with the central star since otherwise the color temperature of this 
emission should fall off rapidly with distance from the illuminating star; instead, 
she proposed that this emission is emitted by ultrasmall grains, 10 A in radius, 
which undergo large excursions in temperature due to stochastic heating by 
single stellar photons. 

The presence of a population of ultrasmall grains in the diffuse ISM was 
explicitly indicated by the 12 jim and 25 //m "cirrus" emission detected by the 
Infrared Astronomical Satellite (IRAS) (Boulanger k. Perault 1988), which is far 
in excess (by several orders of magnitude) of what would be expected from large 
grains of 15-25 K in thermal equilibrium with the general interstellar radiation 
field. Subsequent measurements by the Diffuse Infrared Background Experiment 
(DIRBE) instrument on the Cosmic Background Explorer (COBE) satellite con- 
firmed this and detected additional broadband emission at 3.5 /Ltm and 4.9 /Ltm 
(Arendt et al. 1998). 

More recently, spectrometers aboard the Infrared Telescope in Space (IRTS) 
(Onaka et al. 1996; Tanaka et al. 1996) and the Infrared Space Observatory (ISO) 
(Mattila et al. 1996) have shown that the diffuse ISM radiates strongly in emis- 
sion features at 3.3, 6.2, 7.7, 8.6, and 11.3 ^m. These emission features, first seen 
in the spectra of the planetary nebulae NGC 7027 and BDH-30°3639 (Gillett, 
Forrest, k, Merrill 1973), have been observed in a wide range of astronomical en- 
vironments including planetary nebulae, protoplanetary nebulae, reflection neb- 
ulae, HII regions, circumstellar envelopes, and external galaxies (see Ticlcns ct 
al. 2000 for a review). Often referred to as "unidentified infrared" (UIR) bands, 
these emission features arc now usually attributed to frcc-flying PAH molecules 
which are vibrationally excited upon absorption of a single UV/visible photon 
(Leger & Puget 1984; Allamandola, Tielens, & Barker 1985; AUamandola, Hud- 
gins, & Sandford 1999; Draine k Li 2001; Li k Draine 2001b) although other car- 
riers have also been proposed such as carbonaceous grains with a partly aromatic 
character (e.g., HAG [Duley k Wilhams 1981, Duley, Jones, k Wilhams 1990], 
quenched carbonaceous composite [QCC; Sakata et al. 1992], coal [Papoular et 
al. 1996], fullerenes [Geo; Webster 1993]), and surface-graphitized nanodiamonds 
(Jones k d'Hendecourt 2000). 

Since the late 1970s, various modern interstellar dust models have been 
developed (see Li k Greenberg 2003, Draine 2004, Dwek et al. 2004 for recent 
reviews). In general, these models fall into three broad categories: 

• The silicate-graphite model - this model consists of two (physically) sep- 
arate dust components each with a power-law size distribution dn{a)/da ~ 



Interstellar Nanoparticles 



5 



a in the size range 50 A ^ a ^ 0.25 (Mathis, Rumpl, Sz Nordsieck 
1977 [hereafter MRN]; Draine & Lee 1984). Modifications to this model 
were later made by Draine & Anderson (1985), Weiland et al. (1986), 
Sorrell (1990), Siebenmorgen & Kriigel (1992), Rowan-Robinson (1992), 
Kim, Martin, & Hendry (1994), Dwek et al. (1997), Clayton et al. (2003), 
and Zubko, Dwek, & Arendt (2003) by including new dust components 
(e.g., amorphous carbon, carbonaceous organic refractory, and PAHs) and 
adjusting dust sizes (e.g., deriving dust size distributions using the "Max- 
imum Entropy Method" or the "Method of Regularization" rather than 
presuming a certain functional form). 

Recent developments were made by Draine and his coworkers (Li &; Draine 
2001b, 2002a, 2002b; Weingartner & Draine 2001a) who have extended 
the silicate-graphite grain model to explicitly include a PAH component 
as the small-size end of the carbonaceous grain population. The PAH 
component, containing ~ 45 x 10~^ of C relative to H, is represented by 
a log-normal size distribution dn[a)/dlna ~ exp |— [In (a/ag)]^ / (2(T^) | 

with ao « 3.5 A and a « 0.4 for a;^3.5A (see Li & Draine 2001b). Note 
that with the PAH component added, the size distribution of the larger 
grains is changed from the simple MRN power-law (see Weingartner &i 
Draine 2001a; Zubko et al. 2003). 

The silicate core-carbonaceous mantle model - originally proposed by 
Greenberg (1978) in the context of a complex cyclic evolutionary scenario 
(see Greenberg & Li 1999), this model consists of larger silicate grains 
coated by a layer of carbonaceous organic refractory material, produced 
by UV photolysis of ice mixtures which attempts to simulate the physical 
and chemical processes occurring in interstellar space. The most recent 
development of this model was that of Li & Greenberg (1997), who mod- 
eled the core-mantle grains as 2:1 (the ratio of the length to the diameter) 
finite cylinders (to account for the interstellar polarization) with a Gaus- 
sian size distribution for the mantle dn{a) /da ~ exp —5 (a — Oc)^ /af and 
a single silicate core radius ac~0.07 /um and a "cut-off" size ai~0.066 fim. 
In addition, a PAH component and a population of small graphitic grains 
are added respectively to account for the far-UV extinction rise plus the 
"UIR" emission bands and the 2175 A extinction hump. 

Again, modifications to this model were also made by considering dif- 
ferent coating materials (e.g., amorphous carbon, HAG), including new 
dust type (e.g., iron, small bare sihcates), and varying dust size distribu- 
tions (Chlewicki & Laureijs 1988; Duley, Jones, & Williams 1989; Desert, 
Boulanger, & Puget 1990; Li & Greenberg 1998; Zubko 1999a). The Duley 
et al. (1990) model exclusively consists of silicate core-HAC mantle grains 
with a bimodal size distributions: (1) ^ 43% of the total dust mass in 
grains ^100 A which produce the far-UV extinction rise and the 2175 A 
hump through an electronic transition of the 0H~ ions in low-coordination 
sites on or within silicate grains, and (2) the remaining dust mass in grains 
0.05^a^0.25 /im with a power-law dn{a)/da ~ a~^'^ size distribution. The 
"UIR" bands were attributed to the thermally isolated aromatic "islands" 
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of HAC material - Duley et al. (1989) postulated that an absorbed stellar 
photon might remain localized in a single aromatic unit of the grain long 
enough for the unit to cool by IR vibrational relaxation rather than by 
transferring the photon energy to the phonon spectrum of the grain. 

• The; composite grain model - realizing that grain shattering due to grain- 
grain collisions and subsequent reassembly through agglomeration of grain 
fragments may be important in the ISM, Mathis &l Whiffen (1989) modeled 
the interstellar grains as composite collections of small silicates, vacuum 
(~ 80% in volume), and carbon of various kinds (amorphous carbon, HAC, 
organic refractories) with a power-law size distribution dn{a) / da^ a~^'^ in 
the size range 0.03 /xm ;^ a ^ 0.9 /xm. In addition, a separate small graphite 
component containing « 59x 10~^ of C is needed to account for the 2175 A 
extinction hiimp. Wright (1987) argued that a fractal structure would be 
expected for interstellar grains formed through the coagulation of small 
grain fragments created from grain disruption caused by supernova shock 
waves. In view that the relative abundances of refractory elements in the 
ISM may be as low as 65% of solar (Snow & Witt 1996), Mathis (1996) 
revised the composite model in order to satisfy the tighter abundance con- 
straints. To optimize the use of the heavy elements, Mathis (1996) derived 
a vacuum fraction of ~45%. The "UIR" emission, which remained unac- 
counted for in the models of Mathis & Whiffen (1989) and Mathis (1996), 
was taken into account in Zubko et al. (2003) by including a population 
of PAH molecules. 

It is fair to say that, as can be seen from the above summarizing description 
of the key contemporary dust models, a consensus on the interstellar dust com- 
positions and sizes is now approaching among various grain models, although 
the debate in details is still going on and is not expected to disappear in the 
near future (except that the IR spectropolarimetric observation at 3.4 //m and 
9.7 ^m may allow a direct test of the core-mantle model for interstellar dust; see 
Li & Greenberg 2002) — 

• Regarding the grain chemical composition, the most generally accepted 
view is that interstellar grains consist of amorphous silicates and some 
form of carbonaceous materials; the former is inferred from the 9.7 //m Si- 
O stretching mode and 18 fim. 0-Si-O bending mode absorption features 
in interstellar regions as well as the fact that the cosmically abundant 
heavy elements such as Si, Fe, Mg are highly depleted; the latter is mainly 
inferred from the 2175 A extinction hump (and the ubiquitous 3.4 ^m C-H 
stretching vibrational band) and the fact that silicates alone are not able 
to provide enough extinction. For the carbonaceous component, a wide 
range of dust materials have been suggested including amorphous carbon, 
coal. Ceo, diamond, graphite, HAC, PAHs, organic refractory, and QCC 
(see Pendleton 2004 for a review of the carbonaceous component). 

• It is also generally accepted that, through the analysis of the wavelength- 
dependent interstellar extinction and polarization curves as well as the 

near and mid IR emission, interstellar grain sizes may be separated into 
two domains - (1) the "large" grain component (with radii a>0.025/Ltm; 
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including the "classical" grains [with a^O.l /xm]) which is primarily respon- 
sible for the extinction, polarization and scattering at visible wavelengths 
and the IR emission at A;^60 ^m; and (2) the "very small grain" compo- 
nent (with o<0.025^m) which contributes importantly to the extinction 
in the vacuum- UV and emit strongly in the near and mid IR at A^60 jim. 
(see Figs. 8,16 of Li & Draine 2001b) when transiently heated to high tem- 
peratures during quantized absorption events. 

While the size distribution for the "classical grain" component is relatively 

well constrained by fitting the observed interstellar extinction curve for an 
assumed dust composition, our knowledge of the size distribution dn/da 
for the "very small grain" component is better constrained by the interstel- 
lar near and mid IR emission due to the fact that, for A^^O.l /im, these very 
small grains are in the Rayleigh limit (x = 27ra/A-Cl) and their extinc- 
tion cross sections Cext(O) A) per unit volume V are independent of size, so 
that the observed UV/far-UV extinction curve A{X)/Nii (mag cm~^ H~-^) 
only constrains the total volume Vtot of this grain component: ^(A)/A^h = 
1.086 iVjjVCcxt (a, A) (dn/da) da = 1.086 iVn^^tot (Ccxt/^), where A^h is 
the hydrogen column density. In contrast, the near and mid IR inten- 
sity Ix is sensitive to the grain heat capacity (cx a^) which determines the 
maximum temperature to which the grain can reach when illuminated by 
a radiation field: h = I Cabs(a, A) (dn/da) daj Bx{T[a]) (dP/dT) dT 
where Cabs (O) A) is absorption cross section for a grain of radius a 
at wavelength A (for grains in the Rayleigh limit, Cabs^Ccxt), is the 
Planck function, dP/dT is the dust temperature distribution function. Ev- 
idently, dP/dT is a sensitive function of grain size a (see Draine &; Li 2001, 
Li & Draine 2001b; and §3, Fig. 6 of this paper). 

In the following sections of this review, wc will focus on the ultrasmall grain 
component (a^25nm), with particular emphasis on the photophysical processes 
including the stochastic heating and the vibrational excitation of PAH molecules, 
and the excitation of the photoluminescence of silicon nanoparticles. In astro- 
physical literature, one frequently encounters terms like "ultrasmall grains", 
"very small grains", "large molecules", "tiny grains", and "nanoparticles"; to 
avoid confusion, we emphasize here that they are synonymous. In the following 
we will use the term "nanoparticles" - by "nanoparticles" we mean grains of a 
few angstroms to a few tens of nanometers in radius. 

In §2 we summarize the direct and indirect evidence for the existence of 
nanoparticles in the ISM. The vibrational excitation of nanoparticles is detailed 
in §3. In §4 we discuss the physics regarding the excitation of photoluminescence. 
In §5 we present an overview of the nanoparticle species known or proposed to 
exist in interstellar space. Concluding remarks are given in §6. 

2. Nanoparticles in Interstellar Space 

The discovery of presolar nanodiamonds (see §5.4) and TiC nanocrystals (see 
§5.5) in primitive meteorites implies that there must exist such nano-sized species 
in interstellar space. In addition, as already mentioned in §1, the presence of 
nanoparticles in the ISM is clearly indicated by — 
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(1) The ubiquitous distinctive set of "UIR" emission bands at 3.3, 6.2, 7.7, 
8.6, and 11.3 /im. This emission, accounting for ~ 20% of the total power 
radiated by dust, is closely explained by transiently heated PAHs with an 
abundance C/H?a45 parts per million (ppm) and a log-normal size distri- 
bution peaking at ~ 6 A (with ~ 100 carbon atoms; see Li & Draine 2001b 
for details). 

(2) The mid-IR emission at A ^ 60 jim.. This emission, first discovered by 

the IRAS broadband photometry at 12 and 2h pLva and later confirmed 
by the COBE-DIRBE observations, cannot be explained by large grains 
(a ^250 A) heated by the interstellar radiation field to equilibrium temper- 
atures 15-25 K since the predicted emission intensities from large grains 
are smaller than the IRAS 12 and 25 fira emission intensities by several or- 
ders of magnitude. In the diffuse ISM, the emission at A^60 /xm accounts 
for ~35% of the total power radiated by dust. It is well recognized that 
this emission arises from nanoparticles (a ^250 A) stochastically heated by 
single UV/visible photons to temperatures significantly higher than their 
time-averaged temperatures (see Draine & Li 2001 and references therein). 
Even at 60 /um, the nanoparticle component contributes ^ 70% of the total 
emission power of the diffuse ISM detected by the COBE-DIRBE photome- 
ters (see Fig. 8 of Li & Draine 2001b). 

(3) To a lesser degree, the far-UV extinction rise. The far-UV interstellar 
extinction continues to rise up to A = 0.1 //m, and there does not appear to 
be any evidence of saturation even at this wavelength (see Whittet 2003). 
Since it is generally true that a grain absorbs and scatters light most effec- 
tively at wavelengths comparable to its size A ^ 27ra (Kriigel 2003) , we can 
therefore conclude that there must be appreciable numbers of interstellar 
grains with a ^ 0.1 /im/27r w 16 nm (e.g., in the size distribution derived by 
Weingartner Sz Draine [2001a], grains smaller than 2nm provides 80% 
of the total surface area, although they contain only ~ 6% of the total dust 
mass). However, as remarked earlier (see §1), the far-UV extinction curve 
is not able to tell us the details of the size distribution of the nanoparticle 
component. 

(4) The "anomalous" Galactic foreground microwave emission in the 10 — 
100 GHz region. This emission, discovered unexpectedly in recent exper- 
iments to study the angular structure in the cosmic background radia- 
tion, was found to be positively correlated with interstellar dust, as traced 
by the 100 nm IRAS map or the 140 /xm COBE-DIRBE map (see Draine 
1999 and references therein). However, the dust thermal emission at mi- 
crowave frequencies extrapolated from the 100-3000 //m far-IR emission 
radiated by large grains in thermal equilibrium with the interstellar radi- 
ation field falls far below the observed microwave emission (Draine 1999). 
This "anomalous" emission also significantly exceeds the free-free emis- 
sion from interstellar plasma (Draine & Lazarian 1998a). As described by 
Draine and his coworkers (Draine &: Lazarian 1998b; Draine &: Li 2003), a 
number of physical processes, including collisions with neutral atoms and 
ions, "plasma drag" (due to interaction of the electric dipole moment of 
the grain with the electric field produced by passing ions) , and absorption 
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and emission of photons, can drive nanoparticles to rapidly rotate, with 
rotation rates reaching tens of GHz. The rotational electric dipole emission 
from these spinning nanoparticles, the very same grain component required 
to account for the "UIR" emission and the IRAS 12 and 25 fim emission, 
was shown to be capable of accounting for the "anomalous" microwave 
emission (Draine & Lazarian 1998a,b; Draine 1999). 

• (5) The photoelectric heating of the diffuse ISM. Grains are long thought 
to be an important energy source for the interstellar gas through ejection 
of photoclcctrons since (a) photons with energies below the ionization po- 
tential of H (~ 13.6 eV) do not couple directly to the gas; and (b) other 
heating sources such as cosmic rays, magnetic fields, and turbulence are 
unimportant as a global heating source for the diffuse ISM. 

The photoelectric heating starts from the absorption of a far-UV photon 
by a dust grain, followed by ejection of an electron which then collisionally 
heats the interstellar gas by transferring (to the gas) the excess energy 
left over after overcoming the work function (the binding energy of the 
electron to the grain) and the electrostatic potential of the grain (if it is 
charged). 

In the diffuse ISM, nanoparticles (and in particular, angstrom-sized PAH 
molecules) are much more efficient in heating the gas than large grains 
(see Tielens &; Peeters 2002 and references therein) since (a) the mean free 
path of an electron in a solid is just ~ 10 A and therefore photoelectrons 
created inside a large grain rarely reach the grain surface; and (b) the total 
far-UV absorption is dominated by the nanoparticle component. Recent 
studies show that grains smaller than 10 nm are responsible for ^ 96% of 
the total photoelectric heating of the gas, with half of this provided by 
grains smaller than 15 A (Bakes &; Tielens 1994; Weingartner &; Draine 
2001b). 

• (6) The Extended Red Emission (ERE) ? This emission, first detected in 
the Red Rectangle (Schmidt, Cohen, & Margon 1980), is characterized 
by a broad, featureless band between '^5400 A and 9500 A, with a width 
600 A ^ FWHM <^ 1000 A and a peak of maximum emission at 6100 A ^ 
Ap ^ 8200 A, depending on the physical conditions of the environment 
where the ERE is produced. The ERE has been seen in a wide variety 
of dusty environments: the diffuse ISM of our Galaxy, reflection nebulae, 
planetary nebulae, HII regions, and other galaxies (see Witt &; Vijh 2004 
for a review). The ERE is generally attributed to photoluminescence (PL) 
by some component of interstellar dust, powered by UV/visible photons. 

The observational evidence shows that (a) the ERE carriers must have 
a photon conversion efficiency ?7pl (the number ratio of PL photons to 
exciting photons) substantially larger than 10% as estimated from the cor- 
relation of ERE intensity with HI column density at high Galactic latitudes 
(Gordon, Witt, &: Friedmann 1998), and (b) the carriers can be easily mod- 
ified or destroyed by intense UV radiation (Witt 2000) . This suggests that 
the ERE carriers arc very likely in the nanometer size range because (a) 
in general, nanoparticles are expected to luminesce efficiently through the 
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recombination of the electron-hole pair created upon absorption of an en- 
ergetic photon, since in such small systems the excited electron is spatially 
confined and the radiationless transitions that are facilitated by Auger and 
defect related recombination axe reduced (see §4) ; and (b) small nanoparti- 
cles may be photolytically more unstable and/or more readily photoionized 
in regions where the radiation intensity exceeds certain levels of intensity 
and hardness, and thus resulting in both a decrease in the ERE intensity 
and a rcdshift of the ERE peak wavelength, since (i) photoionization would 
quench the luminescence of nanoparticles, and (ii) the smaller grains would 
be selectively removed due to size-dependent photofragmentation (Smith 
& Witt 2002). 

Proposed ERE carriers include (a) "classic" submicron-sized carbonaceous 
materials: HAC (Duley 1985; Witt & Schild 1988), QCC (Sakata et al. 
1992), and coal (Papoular et al. 1996); (b) nanometer-sized carbon-based 
materials: PAHs (d'Hendecourt et al. 1986), and carbon nanoparticles 
(Seahra & Duley 1999); (c) nanometer-sized silicon-based materials: crys- 
talline silicon nanoparticles (Witt et al. 1998; Ledoux et al. 1998, 2001; 
Smith &; Witt 2002); and (d) particle-bombarded silicate grains (Koike et 
al. 2002). 

The carbon-based models appear to be ruled out: (a) submicron-sized car- 
bon materials appear to be unable to simultaneously match the observed 
ERE spectra and the required PL efficiency (Witt et al. 1998); (b) although 
high photoluminescence efficiencies can be obtained by PAHs, the lack of 
spatial correlation between the ERE and the PAH IR emission bands in 
some regions (see Li &: Draine 2002a and references therein), together with 
nondetection of ERE emission in reflection nebulae illuminated by stars 
with effective temperatures T^,<7000K (Darbon, Perrin, & Sivan 1999), 
whereas PAHs emission bands have been seen in such regions (e.g., see 
Uchida, Sellgren, & Werner 1998) and are expected for the PAH emission 
model (Li k. Draine 2002b), seems to argue against PAHs as ERE carriers; 
(c) the carbon nano-clustcr hypothesis put forward by Seahra &: Duley 
(1999) appears to be invalid as indicated by nondetection in NGC 7023 of 
the predicted 1 /xm ERE peak (Gordon et al. 2000), although they argued 
that these carbon nanoparticles with mixed sp^ /sp^ bonding would be able 
to meet both the ERE profile and the PL efficiency requirements. 

The silicon nanoparticle (SNP) model, originally proposed by Witt et al. 
(1998) and Ledoux et al. (1998), seems promising. Experimental data show 
that SNPs provide so far the best match to the observed ERE spectra 
and to the quantum efficiency requirement. However, this model also has 
difficulty: wc calculated the thermal emission expected from such particles, 
both in a reflection nebula such as NGC 2023 and in the diffuse ISM; we 
found that SNPs would produce a strong emission feature at 20 /7,m which 
is not seen in the observational spectra; therefore we concluded that if the 
ERE is due to SNPs, they must be either in clusters or attached to larger 
grains (see Li & Draine 2002a for details). 
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3. Excitation of Vibrational Transitions 

For isolated, frcc-flying nanoparticles or large molecules in many astrophysi- 
cal environments, photon absorption is often the dominant excitation process. 
Prior to the absorption of an energetic photon, a molecule stays at the ground 
electronic state, which is (the lowest singlet state) for neutrals and Dq (the 
lowest doublet state) for ions, containing very little vibrational energy. After 
photoabsorption occurs, the molecule is excited to an upper electronic state Sn 
for neutrals or Dn for ions (n>l). The electronically excited molecule has 3 
major competing decay channels to relax its energy: radiation, photoionization, 
and photodissociation. In this section we will focus on the radiation process and 
use PAHs as an example. 

3.1. An Overview of the Photoexcitation and Emission Processes 

Once the molecule is electronically excited, both radiative and radiationless 
transitions will occur between its vibrational modes, electronic and vibrational 
states (see Birks 1970 for detailed discussions on neutral molecules): 

• the molecule will undergo internal vibrational redistribution (IVR) , an iso- 
energetic radiationless process between vibrational modes, which rapidly 
(~ 10~^^-10~^'^ s) spreads the absorbed energy among its vibrational de- 
grees of freedom; 

• the molecule will undergo internal conversion (IC), a radiationless transi- 
tion between electronic states of the same multiplicity, which very quickly 
(~ 10~^^-10~^ s) transfers the electronic energy to a highly vibrationally 
excited state of the lower lying electronic state (i.e., Si^Sj for neutrals, 

and Di^Dj for ions where i>j)', 

• the molecule will also undergo intersystem crossing (ISC), a rapid radia- 
tionless transition between electronic states of different multiplicity (i.e., 
Si^Tj for neutrals, and Di^Qj for ions where i>j); 

• the molecule will undergo fluorescence , a radiative electronic transition 
(~ 10~^ s) between states of the same multiplicity, which results in the 
emission of a visible photon (i.e., Si^Sj for neutrals, and Di^Dj for ions 
where i>j; for many molecules, this transition mainly takes place between 
the first excited electronic state or Di and the ground state or Dq, 
as a consequence of rapid internal conversion which results in a practically 
complete conversion of the initial excitation energy to vibrationally excited 
levels of the Si or Di and 5*0 or Do electronic states); 

• the molecule will also undergo phosphorescence with emission of a visible 
photon, a very slow (~ a few seconds) radiative electronic transition be- 
tween states of different multiplicity (i.e., Ti^Sj for neutrals where i>j; 
for ions, this transition [i.e. Qi—^Dj] is not important since internal conver- 
sion dominates over intersystem crossing so that the molecule ends up in 
highly vibrationally excited levels of the ground electronic state Do [Leach 
1987]); 
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• finally, the molecule will undergo IR emission , a radiative vibrational tran- 
sition 0.1 s) between higher and lower vibrational states of the same elec- 
tronic state. This transition mainly takes place in the ground electronic 
state Sq or Dq because of rapid internal conversion and electronic fluores- 
cence which ultimately drive the molecule down to the ground electronic 
state with a high vibrational energy, except for some neutral molecules 
(e.g. chrysene C18H12) a large fraction of this transition occurs in the first 
excited triplet state Ti as a result of rapid intcrsystcm crossing. Note that 
in Figure 1 we only plot the Av=l vibrational transitions (where v is the 
vibrational quantum number) given by the harmonic oscillator selection 
rule; but overtone transitions with Av=2,3,... are also allowed when ac- 
count is taken of anharmonicity. We should stress that in astrophysical 
literature, vibrational transitions are sometimes also called "vibrational 
fluorescence" (e.g. see Allamandola et al. 1989); in this case, we shall call 
the fluorescence process described above "electronic fluorescence" . 

In addition, for highly isolated molecules, there could exist another two ra- 
diative energy decay channels - the inverse fluorescence (Leach 1987) and the 
recurrent fluorescence (Leger, Boissel, & d'Hendecourt 1988). The former re- 
sults from a transition from a high vibrational level of a lower electronic state to 
a higher electronic state. The latter, also known as the Poincare fluorescence , 
was postulated by Leger et al. (1988) as resulting from an inverse electronic 
conversion (i.e., the partial conversion of the vibrational energy of the ground 
electronic state into electronic excitation), followed by emitting a visible photon 
just like the ordinary (electronic) fluorescence process. However, the Poincare 
fluorescence may have a quantum yield larger than one (i.e., several fluorescence 
photons can be emitted during one photoabsorption event since the molecule 
can oscillate many times between the electronic ground state and the excited 
state provided the absorbed energy is high enough; see Leger et al. 1988). 

Figures 1 and 2 give a schematic overview of the above-described main pro- 
cesses following the absorption of an energetic photon respectively for a neutral 
molecule and a molecular cation. 

The stochastic heating of ultrasmall grains has been extensively studied 
in literature (see Draine &: Li 2001 and references therein). All studies prior 
to the identification of PAH molecules as the carrier of the "UIR" bands used 
the "thermal" approach (i.e., the vibrationally excited grain is described as 
a system having an internal temperature; see §3.3). The issue of "thermal" 
versus "non-thermal" arose when Leger and his coworkers used the "thermal" 
approach to model the vibrational excitation of PAHs to explain the "UIR" 
bands (Leger k. Puget 1984; Leger, d'Hendecourt, k, Defourneau 1989), while 
Allamandola, Tielens, & Barker (1985, 1989) took a statistical approach and 
questioned the validity of the "thermal" approximation. Barker & Cherchneff 
(1989), d'Hendecourt et al. (1989), Schutte, Tielens, & Allamandola (1993), and 
Cook & Saykally (1998) found that the thermal approximation was valid for 
computing thermal emission spectra. Below we will discuss the recent efforts 
we took to model the PAH excitation and de-excitation processes (Draine & 
Li 2001) using both the "exact-statistical" approach (§3.2) and the "thermal- 
approximation" approach (§3.3). 
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Figure 1. Step ladder model of a large neutral molecule. The 
schematic energy levels involve an electronic part (singlets 5o, Si, S2 
... and triplets Ti, T2 ...) and a vibrational part (vi, V2, V3...). Upon 
absorption of an UV/visible photon, the molecule makes an electronic 
transition from the ground state Sq to an upper state Sn, followed by 
(1) radiationless internal vibrational redistribution, (2) radiationless 
electronic transitions (internal conversion [IC] and intersystem cross- 
ing [ISC]), (3) radiative electronic transitions (fluorescence and phos- 
phorescence), (4) radiative inverse electronic transitions (inverse fluo- 
rescence, and Poincare fluorescence resulted from the inverse internal 
conversion [IIC]), and (5) IR vibrational transitions. 
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Vibrational Excitation/De-excit. of Ions 




Figure 2. Same as Figure 1 but for a large molecular cation (with 
one unpaired electron). The schematic energy levels involve an elec- 
tronic part (doublets Dq^ Di, D2 ... and quartets Qq, Qi, Q2 ■■■) and a 
vibrational part (vi, V2, V3...). In comparison with neutral molecules, 
radiative electronic transitions (fluorescence and phosphorescence) for 
ions are not important since internal conversion quickly transfers al- 
most all of the initial excitation energy to the vibrationally excited 
ground electronic state Dq (see Leach 1987; Allamandola et al. 1989). 
Therefore, IR emission becomes almost the only deactivation route. 
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3.2. "Exact-Statistical" Approach 

As discussed in §3.1, in most cases soon after the photoabsorption an isolated 
nanoparticle (or large molecule) converts almost all of the initial photoexcitation 
energy to vibrational energy of the highly vibrationally excited ground electronic 
state, and hence for both neutrals and ions, IR emission is always the dominant 
deactivation process. Therefore, it is reasonable to model the stochastic heating 
of a nanoparticle in terms of pure vibrational transitions. 

Ideally, if both the vibrational energy levels and the level-to-level transition 
probabilities were known, we could (at least in principle) solve for the statistical 
steady-state populations of the different energy levels of grains illuminated by 
a known radiation field. However, this level of detailed information is generally 
unavailable, for even the smallest and simplest PAH molecules. 

Draine Sz Li (2001) developed an "exact-statistical" theory for modeling the 
photoexcitation and emission processes of nanoparticles. In this theory, the state 
of the grain is characterized by its vibrational energy E. Since there are too many 
vibrational energy levels to consider individually, they are grouped into (M -|- 1) 
"bins" j = 0, ...,M, where the j-th bin is [-Ej^min, ^j,max)) with representative 
energy Ej = {Ej^„im + £'j,max)/2, and width AEj = (£'j,max - ^i,min) (see Fig. 3 
for illustration). Let Pj be the probability of finding the grain in bin j. The 
probability vector Pj evolves according to dPjdt = '^ijPj ^ ^j^i '^ ji-^i 

i = 0, 1, M, where the transition matrix element is the probability per unit 
time for a grain in bin j to make a transition to one of the levels in bin i. If we 
define the diagonal elements of T to be Tjj = — Y^j^i T^j, then under the steady 
state condition (i.e. dPi/dt = for i = 0, 1, M) the state probability evolution 
equation becomes Yj^^TijPj = for i = 0, ...,M. Combining this with the 
normalization condition J2jLoPj = 1, we obtain a set of M linear equations for 
the first M elements of Pj: J2f=o^ {T^ij - T^m) Pj = -TiM for i = 0, M - 1, 
which we solve using the bi-conjugate gradient (BiCG) method. The remaining 
undetermined element Pm is obtained by Pm = — (Tmm)^^ ^j=o^ '^MjPj- 

For a given starlight energy density ue, the state-to-state transition rates 
Tji for transitions i— >j can be determined from photon absorptions and pho- 
ton emissions. The rate for upward transitions l—^u is just the absorption 
rate of photons with such an energy that they excite the grain from bin / 
just to bin u. If the bin width is sufficiently small (i.e., if max [A£^/, AE^^j] <C 
[Eu — El]), the l^u excitation rate is simply Ty^i « Cabs(-^) cue^Eu/ {Eu — Ei) 
for u<M, where CahsiE) is the grain absorption cross section at wavelength 
X=hc/E {h is the Planck function and c is the speed of light); and Tmi ~ 
Cabs(^) cueA£;m/ {Em - Ei) + J^^_E^ dEC^UE)cuE/ {Em - Ei), where the 
integral takes energy absorbed in transitions to levels beyond the highest bin and 
allocates it to the highest bin (M). For the special case of transitions u—l^u we 

include "intrabin" absorptions: /o^^"~^ dE ^1 - -^^—^ Cabs{E) cue/ {Eu — Eu-ij- 
Correction for finite bin width, which is important when the treatment is ap- 
plied to grains with radii a^50A, has been made by Draine & Li (2001) by 
introducing a Gui{E) factor (see Eqs. [15-25] of Draine & Li 2001). 

The rates for downward transitions -u— >Z can be determined from a de- 
tailed balance analysis of the Einstein A coefficient. Similarly, if the bin width 
is sufficiently small, the u^l de-excitation rate can be approximated as T/^ f» 
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Figure 3. Schematic diagram of the vibrational excitation and re- 
laxation processes in nanoparticles (or large molecules). The vibra- 
tional energy levels are grouped into {M + 1) "bins" j = 0,...,M, 
where the j-th bin is [i^'j,mini -E'j^max), with representative energy Ej 



Ej,mm)- The ground 



{Ej,rain + -E'j,max)/2, and width I\Ej 

state is set at -Ecmin = -E'o.max = £"0 = (see Appendix B of Draine & 
Li 2001 for the procedure for specifying the bins). The wiggly arrow 
represents an event of photon absorption {left) or emission [right). The 
small arrows within the bin u and bin I represent the "intrabin" transi- 
tions. With the bin-to-bin transition rates determined, we can solve for 
the statistical steady-state populations of the different energy states of 
grains illuminated by a known radiation field, and then calculate the 
resulting IR emission spectra. 
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W^;^^-^^^<^abs(£^) '^ + ^wUe for Z<u-1, where the UE-containing term 

is the contribution of stimulated emission, and the degeneracies and gi are 
the numbers of energy states in bins u and I, respectively: gj = iV(£'j^max) ^ 

N{Ej,mhi) ^ (f^) ^Ejj where (4^) is the vibrational density of states at 



^i,minj ^ {lEjj^. ^^i' {dEj 

internal energy Ej, which corresponds to the number of ways of distributing 
this energy between different modes of this grain. Again, we refer the reader to 
Draine & Li (2001) for finite bin width corrections as well as "intrabin" radiation 
consideration (see Eqs. [29-31] of Draine & Li 2001). 

It is seen from the above discussions that we require only Cahs{E), the 
degeneracies gj, and the starlight spectrum ue to completely determine the 
transition matrix Tjj . A molecule containing Na atoms will have = 3 A^a — 6 
distinct vibrational modes (plus 3 translational degrees of freedom and 3 ro- 
tational degrees of freedom). If the molecule is approximated as a set of N„i 
harmonic oscillators, and the frequencies of all normal modes of this molecule 
are known, we can calculate N{E), the number of distinct vibrational states 
with total vibrational energy less than or equal to E, using the Beyer-Swinehart 
algorithm (Beyer &; Swinehart 1973; Stein &; Rabinovitch 1973). So far, the 
frequencies of these normal modes have been computed only for a small number 
of PAHs, with some frequencies determined experimentally, but mode spectra 
are not yet available for most PAHs of interest. 

Since exact densities of states are often unknown for interstellar PAHs, the 
Whitten &; Rabinovitch (1963) approximation, a semi-empirical expression, has 

been extensively used in literature: p{E) = (E + ^ S^)^™"V {Nm - 1)! riil™ hi^i 
where p{E) is the density of states (the number of accessible vibrational states 
per unit energy) at internal energy E, E^ = Xli=L {hi'i/2) is the total zero point 
energy of the molecule, Ui is the vibrational frequency, and 0<^<1 is an empirical 
correction factor. 

In contrast, the Draine &: Li (2001) "exact-statistical" theory does not need 
this approximation; instead, they calculate the "theoretical" mode spectrum 
from the Debye model. A PAH molecule containing Nq C atoms and A^h H 
atoms is treated by Draine &l Li (2001) as having 5 different types of vibration: 
(1) (Ac — 2) out-of-plane C-C modes at X^q^^ = k^@op/hc ^ (16.7//m)^ k, 
600 cm~^ given by a two-dimensional Debye model with a Debye temperature 
©op ~ 950 K, where k-Q is the Boltzmann constant, (2) 2 (A^c ~ 2) in-plane 
C-C modes at \qqi^ = k^Qip/hc « (5.7/im)~^ 1740 cm^^ given by a two- 
dimensional Debye model with a Debye temperature Oip 2500 K, (3) Ah out- 
of-plane C-H bending modes at A^Hop — (ll-3jum)~^ k. 886 cm~^, (4) Ah in- 
plane C-H bending modes at A^h ip = (8.6 jLtm)~^ 1161 cm~^, and (5) Ah C-H 
stretching modes at A^^ ^^.j, = (3.3;um)~-^ ~ 3030 cm~^. The "synthetic" mode 
spectrum for C24H12 is in excellent agreement with the actual mode spectrum of 
coronene (see Fig. 1 of Draine &; Li 2001). Similarly, a silicate grain containing 
Na atoms is treated as having 2 ( A^ — 2) vibrational modes distributed according 
to a two-dimensional Debye model with a Debye temperature = 500 K, and 
{Na — 2) modes described by a three-dimensional Debye model with = 1500 K. 
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Figure 4. The cumulative energy probability distributions for PAHs illuminated 
by the general ISRF computed using the exact- statistical model, the thermal-discrete 
model, and the thermal- continuous model. Note that the lowest energy state [E — 
0), not shown here, has P(E > 0) = 1. Taken from Draine & Li (2001). 

From the "synthetic" model mode spectrum we can obtain the vibrational 
density of states and hence the "degeneracy" gj, the number of distinct quantum 
states included in bin j. We note that the densities of states computed for 
C24H12, using both the actual normal mode spectrum for coronene and the 
model normal mode spectrum for C24H12 (see Fig. 1 of Draine & Li 2001) are 
essentially identical for E/hc ^ 300 cm~^ (see Fig. 4 of Draine &: Li 2001). With 
gj derived from the model mode spectrum, and the j— >-i ij<i) excitation rates 
Tjj calculated from a known radiation field with energy density ue, we can 
determine the z— >j {i>j) de-excitation transition rates Tjj. Solving the steady- 
state state probability evolution equation Ylj^i'^ijPj = S 

j^iTjiPi for i = 

0, 1, M, we are able to obtain the steady-state energy probability distribution 
Pj and calculate the resulting IR emission spectrum. 

In Figure 4 we present the cumulative energy probability distributions for 
selected PAHs excited by the general solar neighbourhood interstellar starlight 
radiation field (ISRF) of Mathis, Mezger & Panagia (1983, hereafter MMP) 
obtained from the "exact-statistical" model. It is seen in Figure 4 that the 
probability of being in the ground state is very large for small grains: for ex- 
ample, for the MMP radiation field, grains with Nq <^ 4000 spend most of their 
time at = 0. The sharp drop at 13.6 eV {E/hc ~ 1.1 x lO^cm"^) is due to 
the radiation field cutoff at 912 A and to the fact that multiphoton events are 
rare. The resulting IR emission spectra are displayed in Figure 5. The sawtooth 
features seen at long wavelengths are due to our treatment of transitions from 
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Figure 5. IR emissivities (per C atom) for selected ionized PAHs in the general 
ISRF calculated using the exact-statistical and thermal-discrete models. Taken from 
Draine & Li (2001). 

the lower excited energy bins to the ground state and first few excited states 
(see §5.1 and Appendix B of Draine & Li 2001). 

3.3. "Thermal" Approach 

Given (1) the Debye model assumption about the vibrational mode spectrum, 
(2) the assumption that the absorption cross section Cabs(-£') depends only on 
the photon energy, and (3) the assumption of rapid internal vibrational redis- 
tribution, the "statistical" approach described in §3.2 is "exact". It is often 
desirable, however, to use an alternative "thermal" approximation which is less 
computationally demanding, since the vibrational density of states p{E) is a 
steeply-increasing function of the vibrational energy E (see Fig. 4 of Draine & 
Li 2001), and for large molecules the number of states is often too large to be 
tractable even when E is small. 

The only difference between the "thermal" approach and the "statistical" 
approach concerns the calculation of the downward transition rates T/^ which, 
in contrast to the "exact-statistical" treatment (see §3.2), uses the notion of 
"grain temperature" : instead of using degeneracies gu and gi , the thermal ap- 
proach replaces gi/gu by exp(fez//fcB'ru)-i ' "^^^^^ ^« is the "temperature" of 
the upper level u; i.e., the thermal approach assumes that the emission from 
a molecule with vibrational energy E and Na atoms can be approximated by 
the average emission of a system containing 3{Na — 2) vibrating harmonic os- 
cillators, each with a fundamental frequency i^j, connected to a heat bath with 
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Figure 6. The energy probability distribution functions for charged carbona- 
ceous grains (a = SA [C60H20], 10 A [C480 H120], 25 A [C7200H1800], 50, 75, 100, 150, 
200, 250, 300 A) illuminated by the general ISRF. The discontinuity in the 5, 10, 
and 25 A curves is due to the change of the estimate for grain vibrational "temper- 
ature" at the 20th vibrational mode (see Draine & Li 2001). For 5, 10, and 25A a 
dot indicates the first excited state, and Pq is the probability of being in the ground 
state. Taken from Li & Draine (2001b). 

temperature T such that the expectation value of the energy of the vibrational 
system would be E{T) = Y^^=i exp(fei/^/fcBT')-i • When the number of modes is 
large, this summation contains many terms. Therefore, we take the Debye model 
discussed in §3.2 for silicate grains and the C-C modes of PAHs to estimate the 
grain "temperature" associated with internal energy E. 

So far, the grains are assumed to undergo "discrete cooling"; i.e., the grains 
make discrete transitions to energy levels / < n by emission of single photons. 
There are substantial computational advantages (Guhathakurta & Draine 1989; 
Draine & Li 2001) if the cooling of the grains is approximated as continuous 
rather than discrete, so that the only downward transition from a level u is to 
the adjacent level u — 1 (i.e. T^^ = for u — I > 1). We refer to this as the 
"thermal continuous" cooling approximation. 

In Figure 4 we also plot the cumulative energy probability distributions for 
selected PAHs excited by the MMP ISRF obtained from the thermal-discrete 
model, and the thermal-continuous model. It is seen that the energy distribu- 
tions P{E) found using the thermal-discrete model and the thermal-continuous 
model are both in good overall agreement with the results of the exact-statistical 
calculation. The IR emission spectra obtained from the thermal-discrete model, 
as shown in Figure 5, are almost identical to those of the exact-statistical model. 
The thermal-continuous cooling model also results in spectra which are very 
close to those computed using the exact-statistical model (not shown here; see 
Figs. 14,15 of Draine & Li 2001). 
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Figure 7. Infrared emission spectra for small carbonaceous grains of various 
sizes heated by the general ISRF, calculated using the full energy distribution func- 
tion P(E) (sohd lines); also shown (broken lines) are spectra computed for grains 
at the "equihbrium" temperature T. Transient heating effects lead to significantly 
more short wavelength emission for a < 200 A. Taken from Li & Draine (2001b). 



Figure 6 shows the energy probabihty distribution functions found for PAHs 
with radii a = 5, 10, 25, 50, 75, 100, 150, 200, 300A illuminated by the general 
ISRF. It is seen that very small grains (a^lOOA) have a very broad P{E), and 
the smallest grains (a<30 A) have an appreciable probability Pq of being found 
in the vibrational ground state E = 0. As the grain size increases, P{E) be- 
comes narrower, so that it can be approximated by a delta function for a>250 A. 
However, for radii as large as a=200 A, grains have energy distribution functions 
which are broad enough that the emission spectrum deviates noticeably from the 
emission spectrum for grains at a single "steady-state" temperature T, as shown 
in Figure 7. For accurate computation of IR emission spectra it is therefore im- 
portant to properly calculate the energy distribution function P{E), including 
grain sizes which are large enough that the average thermal energy content ex- 
ceeds a few eV. 



4. Excitation of Photoluminescence 

Many materials are capable of emitting visible luminescence when subjected 
to some form of excitation such as UV light (photoluminescence), nuclear ra- 
diation such as 7 rays and a and /3 particles (scintillation), mechanical shock 
(triboluminescence), heat (thermoluminescence), chemical reactions (chemilu- 
minescence), and electric fields (electroluminescence). In this section we will 
restrict our attention to photoluminescence, since the interstellar dust lumines- 
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cence, as manifested by the Extended Red Emission ubiquitously seen in inter- 
stellar environments, is believed to be a photon-driven process (Witt & Schild 
1985; Duley 1985; Smith & Witt 2001; Witt & Vijh 2004). 

There are two pre-requisites for luminescence: (1) the luminescent material 
must have a semiconductor structure with a nonzero band gap (e.g. metals 
do not luminesce since they have no band gap); (2) energy must be imparted to 
this material before luminescence can take place. The mechanism of photolumi- 
nescence in semiconductors is schematically illustrated in Figure 8, which plots 
the E-k diagrams for a direct band gap material {left) and an indirect gap ma- 
terial {right), where E and k are respectively the kinetic energy and wave vector 
(or "momentum vector") of the electron or hole {E = /i;^7;,^/2m*, where %=h/2TT 
is the Planck constant h divided by 2tt, and is the electron or hole effective 
mass). The direct and indirect gap materials are distinguished by their relative 
positions of the conduction band minimum and the valence band maximum in 
the Brillouin zone (the volume of k space containing all the values of A: up to 7r/a 
where a is the unit lattice cell dimension). In a direct gap material, both the 
conduction band minimum and the valence band maximum occur at the zone 
center where fc = 0. In an indirect gap material, however, the conduction band 
minimum does not occur at A; = 0, but rather at some other values of k which 
is usually at the zone edge or close to it (see Fox 2001). 

Upon absorption of an UV or visible photon with an energy ^Wexc exceeding 
the band gap E^ (the gap in energy between the valence band and the conduction 
band) of the material, an electron-hole pair is created and the electron (hole) is 
excited to states high up in the conduction (valence) band (see Fig. 8) . 

During a photon absorption process in semiconductors, we must conserve 
both energy and momentum. In a direct band gap material, the conduction 
band minimum and the valence band maximum have the same k values (i.e., 

Tikj = hkf, where kj and kj are respectively the wave vectors of the initial and 
final electron states; this implies that the electron wave vector should not change 
significantly during a photon absorption process) , conservation of momentum is 
guaranteed for the photoexcitation of the electron which only involves a UV 
or visible photon: ^k, -|- Tikphot ^ = Tikf, since kphot; the wave vector of 
the absorbed photon (which is in the order of 27r/A 10"^cm~^), is negligible 
compared to the electron wave vector (which is related to the size of the Brillouin 
zone ir/a ^ 10*^ cm^^, where the unit cell dimension a is in the order of a few 
angstroms). This implies that in a direct band gap material, the electron wave 
vector does not change significantly during a photon absorption process. We 
therefore represent photon absorption processes by vertical arrows in the E-k 
diagrams (see Fig. 8). 

In contrast, for an indirect band gap material of which the conduction 
band minimum and the valence band maximum have different k values (see 
Fig. 8) , conservation of momentum implies that the photon absorption process 
must be assisted by either absorbing or emitting a phonon (a quantum of lattice 
vibration), because the electron wave vector must change significantly in jumping 
from the valence band in state (-Bj,kj) to a state {Ef,hf) in the conduction band, 
and the absorption of a photon alone can not provide the required momentum 
change since |kphot| |kj — k/|. 
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The excited electron and hole will not remain in their initial excited states 
for very long; instead, they will relax very rapidly (~10^^^s) to the lowest 
energy states within their respective bands by emitting phonons. When the 
electron (hole) finally arrives at the bottom (top) of the conduction (valence) 
band, the electron-hole pair can recombine radiatively with the emission of a 
photon (luminescence), or nonradiatively by transferring the electron's energy 
to impurities or defects in the material or dangling bonds at the surface. 

Just like the photon absorption process discussed above, the electron-hole 
recombination in a direct band gap material does not involve any phonons since 
there is no need for momentum change for the electron. In contrast, in an in- 
direct gap material, the excited electron located in the conduction band needs 
to undergo a change in momentum state before it can recombine with a hole 
in the valence band; conservation of momentum demands that the electron-hole 
recombination must be accompanied by the emission of a phonon, since it is not 
possible to make this recombination by the emission of a photon alone. Com- 
pared to the photon absorption process in an indirect gap material for which 
conservation of momentum can be fulfilled by cither absorption or emission of a 
phonon, in the electron-hole recombination process phonon absorption becomes 
negligible, whereas phonon emission becomes the dominant momentum conser- 
vation mediator because (1) the number of phonons available for absorption is 
small and is rapidly decreasing at lower temperatures, whereas the emission of 
phonons by electrons which are already at a high-energy state is very probable; 
and (2) an optical transition assisted by phonon emission occurs at a lower pho- 
ton energy Eg — /ifphon than the gap energy, whereas phonon absorption results 
in a higher photon energy of at least Eg + hvphon, which can be more readily 
re-absorbed by the semiconductor nanoparticle. But we note that the energy 
of a phonon (/it'phon) is just in the order of ~ 0.01 eV, much smaller than the 
energy of the electron-hole recombination luminescence photon. Also because 
prior to the recombination, the electrons and holes respectively accumulate at 
the bottom of the conduction band and the top of the valence band, the energy 
separation between the electrons and the holes approximately equals to the en- 
ergy of the band gap. Therefore, the luminescence emitted by both types of 
semiconductors occurs at an energy close to the band gap Eg. 

The PL efficiency is determined by the competition between radiative and 
nonradiative recombination. For an indirect gap material, the PL process, which 
requires a change in both energy and momentum for the excited electron and 
hence involves both a photon and a phonon, is a second-order process with a long 
radiative lifetime (~ 10~^ — lO^^s), and therefore a relatively small efficiency 
because of the competition with nonradiative combination. In contrast, in a 
direct gap material, the emission of a PL photon does not need the assistance 
of a phonon to conserve momentum. Therefore, the PL process in a direct 
gap material is a first-order process with a much shorter radiative lifetime 

— s) and a much higher PL efficiency in comparison with an indirect 
gap material. 

However, for particles in the nanometer size domain, we would expect sub- 
stantial changes in both the efficiency and the peak energy of the photolumines- 
cence due to the quantum confinement effect. This can be understood in terms 
of the Heisenberg uncertainty principle. Unlike in bulk materials the electrons 
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Figure 8. Schematic band diagrams for the photoluminescence processes in a 
direct gap material (left) and an indirect gap material [right). The shaded states at 
the bottom of the conduction band and the empty states at the top of the valence 
band respectively represent the electrons and holes created by the absorption of 
photons with energy ?itJoxc>i?g. The cascade of transitions within the conduction 
and valence bands represents the rapid thermalization of the excited electrons and 
holes through phonon emission. In a direct gap material (left), the conduction band 
minimum and the valence band maximum occur at the same k values. Both the 
photon absorption and emission (i.e. the electron-hole recombination) processes can 
conserve momentum without the assistance of phonons, since the momentum of the 
absorbed or emitted photon is negligible compared to the momentum of the electron. 
We therefore represent photon absorption and emission processes by vertical arrows 
on E-k diagrams. In an indirect gap material (right), the conduction band minimum 
and the valence band maximum occur at different k values. As a result, to conserve 
momentum, the photon absorption process must involve either absorption (indicated 
by a "+" sign) or emission (indicated by a "— " sign) of a phonon, while the PL 
process requires the emission of a phonon. Since the energy of a phonon (~0.01 eV) 
is much smaller than the energy of the PL photon, for an indirect gap material, the 
peak energy of the PL also roughly reflects its band gap. 
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and holes are free to move within their respective bands in all three directions, 
in nanoparticles the electrons and holes are spatially confined and hence their 
motion is quantized in all three dimensions. The spatial confinement of a parti- 
cle of mass m to a region in a given direction (say, along the x axis) of length 
Ax would introduce an uncertainty in its momentum Ap^ ^ fi/Ax and increase 
its kinetic energy by an amount Eqq ~ (Apx) /2m ~ ?i^/2m(Ax) . A sim- 
ple particle-in-a-box analysis, using the Schrodinger's equation and the effective 
mass approximation, shows that the ground state quantum confinement energy 

would be Eqc ~ ) , where m* = mlm^/ (m* -|- m^) is the reduced ef- 

fective mass of the electron-hole pair (m* and are respectively the effective 
mass of the electron and hole) (Fox 2001). For nano-sized particles, the quan- 
tum confinement effect becomes significant since the confinement energy Eqq 
would be comparable to or greater than their thermal energy E'th ~ ^kBT at the 
temperature range expected for nanoparticles in the diffuse ISM (see Figs. 7,8 
of Draine & Li 2001 and Fig. 3 of Li & Draine 2002a). For example, a silicon 
grain with m* ~ 0.98me and ~ 0.52me (mg is the free electron mass) smaller 
than ~ 8 nm would exhibit quantum effects at a temperature 100 K (which is 
expected for nano-sized silicon dust in the diffuse ISM; see Fig. 3 of Li & Draine 
2002a) with Eqc « 0.83 (a/nm)"^ eV > Eth « 0.013 (T/lOO K) eV. More de- 
tailed studies show that silicon nanocrystals exhibit an a~^-^^ gap-size depen- 
dence: Eg^ Eo + 1.42 (a/nm)~^'^^ eV, where Eq w l.lTeV is the bulk silicon 
band gap (Delerue, Allan, Sz Lannoo 1993). For nanodiamonds, an gap- 
size dependence was derived from the X-ray absorption spectrum measurements 
(Chang et al. 1999; also see Raty et al. 2003): E^ f=i Eq + 0.38(a/nm)~^ eV, 
where Eq ~ 5.47 cV is the bulk diamond band gap. 

It is apparent, therefore, the quantum confinement effect would lead to a 
progressive widening of the band gap of a nano-sized semiconductor as its size is 
reduced, along with a broadening of the electron-hole pair state in momentum 
space (i.e. an increased overlap between the electron and hole wavefunctions) , 
and a decreasing probability for the pair to find a nonradiative recombination 
center, provided that the surface dangling bonds are passivated which would 
otherwise act as traps for the carriers and quench the PL. While the former 
would shift the PL peak to higher energies, the latter two effects would greatly 
enhance the electron-hole radiative recombination probability and result in a 
higher PL efficiency. This is one of the main reasons why silicon nanocrystals 
are proposed recently by Witt and his coworkers (Witt et al. 1998, Smith Sz 
Witt 2001) and by Ledoux et al. (1998, 2001) as the ERE carrier since they are 
capable of luminescing very efficiently (with essentially 100% quantum efficiency) 
in the energy range (~ 1.4-2.4 eV) over which the ERE has been observed in 
astronomical sources, while bulk silicon, an indirect gap semiconductor with a 
band gap of ^ 1.17eV at T = OK, does not luminesce. 

5. An Overview of Interstellar Nanoparticle Species 

5.1. Nano Carbon Grains: Polycyclic Aromatic Hydrocarbon Dust 

As the most abundant interstellar nanoparticle species, nano carbon grains 
(mainly PAHs), containing ~ 15% of the interstellar C abundance (Li & Draine 
2001b), reveal their presence in the ISM by emitting a prominent set of "UIR" 
bands at 3.3, 6.2, 7.7, 8.6, and 11.3/xm (see §1 and §2). Modern research on 



26 



A. Li 



astrophysical PAHs started in the mid 1980s with the seminal studies of Leger 
k, Puget (1984) and Allamandola et al. (1985) who by the first time exphc- 
itly proposed PAH molecules as the "UIR" band carrier. The PAH model is 
now gaining increasing acceptance because of (1) the close resemblance of the 
"UIR" spectra (frequencies and relative intensities) to the vibrational spectra 
of PAH molecules (see Allamandola &; Hudgins 2003 for a recent review); (2) 
the ability of a PAH molecule to emit efficiently in the "UIR" wavelength range 
following single photon heating (Lcgcr & Puget 1984; Allamandola ct al. 1985, 
1989; Drainc & Li 2001; also sec §3); and (3) the success of the PAH model 
in quantitatively reproducing the observed mid-IR spectra of the Milky Way 
diffuse ISM (Li & Draine 2001b), the quiescent molecular cloud SMC Bl^^l in 
the Small Magellanic Cloud (Li & Draine 2002c), and the "UIR" band ratios for 
a wide range of environments ranging from reflection nebulae, HII regions, pho- 
todissociation regions, molecular clouds in the Milky Galaxy to normal galaxies, 
starburst galaxies, and a Seyfert 2 galaxy (Draine &: Li 2001). 

Recently, the PAH model further gained its strength from a close fit to 
the observed "UIR" bands of vdB 133, a UV-poor reflection nebula, which was 
considered as one of the major challenges to identification of the "UIR" bands 
with PAH molecules (see Uchida, Sellgren, k. Werner 1998), since small, neutral 
PAH molecules have little or no absorption at visible wavelengths and there- 
fore require UV photons for excitation. Li &: Draine (2002b) have shown that 
the "astronomical" PAH model, incorporating the experimental result that the 
visual absorption edge shifts to longer wavelength upon ionization and/or as 
the PAH size increases (see Appendix A2 of Li &; Draine 2002b and references 
therein), can closely reproduce the observed IR emission bands of this reflection 
nebula, and is also able to account for the observed dependence of the 12 /xm 
IRAS emission on the effective temperature of the illuminating star (Sellgren, 
Luan, k, Werner 1990). 

Henning k Schnaiter (1998) suggested that nano-sized hydrogcnatcd carbon 
grains may be responsible for the 2175 A extinction hump as well as the 3.4 ^m 
C-H absorption feature. But this seems to be in conflict with the detection of 
the 3.4 /um feature in regions where the 2175 A hump is not seen, e.g., along the 
sightline toward HD 204827 (Valencic et al. 2003) and the Taurus cloud (Whittet 
et al. 2003). 

Kroto et al. (1985) first proposed that Ceo could be present in the ISM with 
a considerable quantity. This molecule and its related species have later been 
proposed as the carriers of the 2175 A extinction hump, the diffuse interstellar 
bands (DIBs), the "UIR" bands, and the ERE (see Webster 1993 and references 
therein). Foing k Ehrenfreund (1994) attributed the two DIBs at 9577 A and 
9632 A to Cgo- However, attempts to search for these molecules in the UV and 
IR were unsuccessful (Snow k Scab 1989; Somerville k Bellis 1989; Moutou et 
al. 1999; Herbig 2000). These molecules were estimated to consume at most 
<0.7ppm carbon (Moutou et al. 1999). Therefore, Cqo is at most a minor 
component of the interstellar dust family. 

5.2. Nano Silicate Grains 

If nano silicate grains are present in the ISM, single-photon heating by starlight 

(see §3) will cause these grains to radiate in the 10 ^m feature (see, e.g., Draine 
k Anderson 1985, Li k Draine 2001a). The absence of a 10 /xm emission feature 
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in IRAS spectra led Desert et al. (1986) to conclude that not more than 1% of 
Si could be in a<15 A silicate grains, and on this basis recent grain models (see 
Li & Draine 2001a and references therein) have excluded silicate nanoparticles 
as a significant grain component. Nondetection of the 10 /^m silicate emission 
feature in the diffuse ISM by ISO (Mattila et al. 1996) and the IRTS (Onaka et 
al. 1996) appeared to confirm this conclusion. 

However, the presence of nano silicate grains can not be ruled out since the 
10 /im silicate emission feature may be hidden by the dominant PAH features. 
Li & Draine (2001a) have placed quantitative upper limits on the abundances 
of both amorphous and crystalline silicate nanoparticles by calculating spectra 
for such tiny grains heated by starlight, and comparing to measurements of the 
IR emission of the diffuse ISM by IRTS (Onaka et al. 1996) and by the DIRBE 
instrument on the COBE satellite. The interstellar extinction curve and the 
10 /um silicate absorption profile are also invoked to provide further constraints. 
It was found that, contrary to previous work, as much as ~10% of interstellar Si 
could be in a^lSA silicate grains without violating observational constraints, 
and not more than ~5% of the Si can be in crystalline silicates (of any size). 

5.3. Nano Silicon Grains 

The presence of a population of silicon nanoparticles containing ^5% of the 
total interstellar dust mass (with Si/H « 6ppm) in the ISM was proposed by 
Witt and his coworkers (Witt et al. 1998, Smith & Witt 2002) and Ledoux et al. 
(1998) to account for the ERE phenomenon. Witt et al. (1998) suggested that 
the formation of interstellar SNPs could occur as a result of the nucleation of 
SiO molecules in oxygen-rich stellar outflows, followed by annealing and phase 
separation into an elemental silicon phase in the core and a passivating mantle 
of Si02. The Si02 mantle is crucial for the SNP model since SNPs luminesce 
efficiently only when their surface dangling Si bonds are passivated. 

The SNP model is recently receiving much attention because (1) as a conse- 
quence of quantum confinement (see §4), the SNP photoluminescence provides 
so far the best match to the observed ERE in both spectral profile (i.e. peak 
position and width) and the required extremely high quantum efficiency: (2) as 
shown by Smith & Witt (2002) in terms of a SNP photoionization and photofrag- 
mentation theory together with the experimentally established fact that pho- 
toionization quenches the PL of SNPs, this model successfully reproduces the 
observed dependences of the ERE intensity, the ERE quantum efficiency, and 
the ERE peak wavelength on the intensity and hardness of the UV/ visible radi- 
ation fields of a wide variety of dusty environments; (3) as already mentioned in 
§2, other candidate materials all fail to simultaneously satisfy the ERE spectral 
characteristics and the high quantum efficiency requirement; it is worth noting 
that HAC, the long-thought ERE candidate, also falls in this category: HAC 
luminesces efficiently only in the near UV (quantum efficiencies as high as 10% 
have been reported; see Rush, Robertson, &; Amaratunga 1996 and references 
therein) when it has a higher degree of hydrogenation (indicating a larger band 
gap), dehydrogenation is expected to reduce the band gap and thus can redshift 
the PL peak into the wavelength range where the ERE is observed, but this 
results in an exponential drop in its PL efficiency by a factor of ~ 10^. 

As already mentioned in §2, although the SNP model seems very attractive, 
it also has some weakness: it has been shown by Li & Draine (2002a) that SNPs 
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with oxide coatings, if they are free-flying in the ISM, would emit strongly in 
the 20 jJLVti 0-Si-O bending band. Existing COBE-DIRBE 25 jim. photometry ap- 
pears to already rule out such high abundances of SNPs. If SNPs are responsible 
for the ERE from the diffuse ISM, they must either be in a;^50A clusters, or 
attached to larger "host" grains. This problem can not be resolved by invoking 
hydrogen passivation of SNPs because H-passivated SNPs luminesce at blue and 
near-UV wavelengths (Zhou, Brus, & Friesner 2003) while blue PL is not ob- 
served under interstellar conditions (Rush & Witt 1975; Vijh, Witt, & Gordon 
2003). Very recently, Witt k Vijh (2004) argued that Fc- or C-passivated SNPs 
could overcome this problem. It would be very helpful to perform detailed calcu- 
lations of the 11.3 ;Um Si-C emission feature of C-passivated SNPs and compare 
with the upcoming SIRTF [Space Infrared Telescope Facility) high resolution 
spectra of the diffuse ISM. Previously, Whittet, Duley, k Martin (1990) de- 
duced that the abundance of Si in SiC dust is no more than ~ 5% of that in 
silicates, through an analysis of the 7.5-13.5 //m absorption spectra of 10 Galac- 
tic Center sources. It would also be very useful to investigate the IR emission 
properties of Fe-passivated SNPs which luminesce even more efficiently than O- 
passivated SNPs in the wavelength range over which the ERE is observed (Mavi 
et al. 2003). 

5.4. Nano Diamond Grains 

Cosmic nanodiamonds were first detected in primitive carbonaceous meteorites 
and identified as presolar in origin based on their isotopic anomalies (Lewis et 
al. 1987), although their presence in the ISM was proposed almost two decades 
earlier by Saslaw & Gaustad (1969) to explain the interstellar UV extinction 
curve. Five years later, Allamandola et al. (1992) attributed the 3.47 /Ltm ab- 
sorption band seen toward a large number of protostars to the tertiary C-H 
stretching mode in diamond-like carbonaceous materials. Jones k d'Hendecourt 
(2000) further suggested that surface-reconstructed (to sp^-bonded carbon) nan- 
odiamonds could be responsible for the "UIR" emission features, the 2175 A 
extinction hump (also see Sandford 1996), and a part of the far-UV extinction 
at A"^ ;> 7 Hm . More recently, circumstellar nanodiamonds were identified in 
the dust disks or envelopes surrounding two Herbig Ae/Be stars HD 97048 and 
Elias 1 and one post-asymptotic giant branch (AGB) star HR 4049, based on 
the 3.43 /um and 3.53 jj,m C-H stretching emission features expected for surface- 
hydrogenated nanodiamonds (Guillois, Ledoux, k Reynaud 1999; van Kerck- 
hoven, Tielens, k Waelkens 2002). 

Presolar meteoritic nanodiamonds were found to have a log-normal size dis- 
tribution with a median radius ~ 1.3 nm (Lewis, Anders, k Draine 1989) and an 
abundance as much as ~0.1% of the total mass in some primitive meteorites, 
more abundant than any other presolar grains by over two orders of magnitude. 
However, it is important to recognize that this abundance may not be represen- 
tative of its original ISM proportion but rather its durability (see Draine 2003a) . 
Indeed, Lewis et al. (1989) found that in interstellar space, as much as 10% of the 
interstellar carbon (~ 36 ppm) could be in the form of nanodiamonds without 
violating the constraints placed by the interstellar extinction curve. However, 
a much more stringent upper limit of '--'O.lppm was derived by Tielens et al. 
(2000) based on nondetection of the characteristic 3.43 //m and 3.53 //m C-H 
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stretching emission features in the ISM. Of course, interstellar nanodiamonds 
could be more abundant if the bulk of them are not hydrogenated. 

5.5. Nano Titanium Carbide Grains 

Presolar titanium carbide (TiC) grains were first identified in primitive mete- 
orites as nano-sized inclusions embedded in micrometer-sized presolar graphite 
grains (Bernatowicz et al. 1996). Very recently, von Helden et al. (2000) proposed 
that TiC nanocrystals could be responsible for the prominent 21 fim emission 
feature detected in over a dozen carbon-rich post-AGB stars. This mysterious 
21 /im feature still remains unidentified since its first detection in the IRAS LRS 
(Low Resolution Spectrometer) spectra of four post-AGB stars (Kwok, Volk, 
&; Hrivnak 1989). While bulk TiC samples show no resonance around 21 /xm 
(Henning Sz Mutschke 2001), laboratory absorption spectra of TiC nanocrystal 
clusters exhibit a strong 21 fim feature, closely resembling the observed 21 /im 
feature both in position, width and in spectral detail (von Helden et al. 2000). 

We have examined the nano-TiC proposal by comparing the maximum 
available with the minimum required TiC dust mass inferred from the nano-TiC 
model (see Li 2003). It is found that if the UV/visible absorption properties of 
TiC nanograins are like their bulk counterparts, the model-required TiC dust 
mass would exceed the maximum available TiC mass by over two orders of mag- 
nitude (also see Chigai et al. 2003; Hony et al. 2003). One may argue that nano 
TiC might have a much higher UV/visible absorptivity so that the available 
TiC mass may be sufficient to account for the observed 21 nm emission feature. 
However, the Kranicrs-Kronig dispersion relations, which relate the wavelength- 
integrated extinction cross sections to the total dust mass (PurccU 1969; Draine 
2003b), would impose a lower bound on the TiC mass. We have shown that 
this Kramers-Kronig lower limit exceeds the maximum available TiC mass by a 
factor of at least ~ 50 (see Li 2003 for details). Therefore, it is unlikely that TiC 
nanoparticles are responsible for the 21 /xm emission feature of post-AGB stars. 

6. Concluding Remarks 

Interstellar nanoparticles are important as emitters of near and mid IR radiation, 
as absorbers of far-UV radiation, as heating agents of the interstellar gas, and 
possibly as luminescing agents of red light. Our understanding of interstellar 
nanoparticles obviously remains incomplete. Besides the long-standing problems 
regarding the carriers of the 2175 A extinction hump and the ERE, there are 
many unanswered questions which will be demanding close attention in the 
future: 

• The origin and evolution of interstellar PAHs are not very clear. Sug- 
gested sources for interstellar PAHs include (1) injection (into the ISM) of 
PAHs formed in carbon star outflows (Latter 1991); (2) shattering of car- 
bonaceous interstellar dust or of photoprocessed interstellar dust organic 
mantles (Greenberg et al. 2000) by grain-grain collisions in interstellar 
shocks (Jones, Ticlens, & Hollcnbach 1996); (3) in-situ formation through 
ion-molecule reactions (Herbst 1991). Although interstellar PAHs contain- 
ing more than ~ 20-30 carbon atoms can survive the UV radiation field 
(Guhathakurta & Draine 1989; Jochims et al. 1994; Allain, Leach, & Sedl- 
mayr 1996; Le Page, Snow, &: Bierbaum 2003), they are still subject to 
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destruction by sputtering in interstellar shock waves (but also see Tielens 
et al. [2000] who argued that destruction of PAHs by sputtering is unim- 
portant, except in extreme environments such as very young supernova 
remnants with ;^200kms~^ shocks, because they couple very well dynam- 
ically to the gas which cools down rapidly), chemical attack by atomic 
oxygen, and coagulation in dense regions (Draine 1994). Clearly, a de- 
tailed study of the evolution of interstellar PAHs would be very valuable. 

The physics and astrophysics of the radiative electronic transitions of in- 
terstellar PAHs are not fully understood. 

Small neutral PAHs with ^40 C atoms are expected to emit near-UV 
and blue photons through fluorescence, phosphorescence, and perhaps also 
through the recurrent Poincare fluorescence. Laboratory studies show that 
the fluorescence quantum yield can be quite high for isolated molecules 
(e.g., fluorescence quantum yields in the range of 10%-45% have been mea- 
sured for gas-phase, collision-free naphthalene CioHg [Reyle & Brechignac 
2000]). However, no PL shortward of A=5000A has been seen in the ISM 
(Rush & Witt 1975; Vijh, Witt, k Gordon 2003). Although PAH ions do 
not exhibit strong ordinary fluorescence or phosphorescence (see the cap- 
tion of Fig. 2), they can undergo Poincare fluorescence which may have a 
quantum yield larger than one. Therefore, it is difficult to explain nonde- 
tection of blue PL in the ISM in terms of PAH photoionization. Moreover, 
for PAHs containing ^40 C atoms in the diffuse ISM, the probability of 
finding them in a nonzero charge state is smaller than ~ 30% (see Fig. 7 
of Li &; Draine 2001b). One may argue that interstellar PAHs are larger 
(which is true, as shown in Li & Draine 2001b, the average [diffuse ISM] 
PAH size is -^6 A, corresponding to 100 C atoms) so that their fluores- 
cence mainly occurs at longer wavelengths, say, in the wavelength range 
over which the ERE is observed. 

However, if interstellar PAHs indeed luminesce in the ERE band with the 
required quantum efficiency (~ 100%; see §2), the "UIR" emission bands 
expected for these luminescing PAHs would be strongly suppressed because 
a considerable fraction of the excitation energy is released in the form of 
PL photons (e.g., in the diffuse ISM with a mean energy Ri2.1eV for 
ERE photons [Szomoru & Guhathakurta 1998] , only ~ 60% of the original 
excitation energy is available as heat for PAHs which have a mean absorbed 
photon energy ~5.2eV [Draine &: Li 2001]; this could be even worse since 
the Poincare fluorescence may have a quantum yield higher than one) . One 
consequence of this would be that even more carbon atoms in the form 
of PAHs would be needed to explain the intensity of the "UIR" bands, 
deteriorating the already tightened carbon budget "crisis" (Snow & Witt 
1996). 

It is unclear how different the optical and thermal properties of nanometer- 
sized materials are compared with their bulk counterparts, although it is 
generally believed that they may be very different. 

— For a small metallic grain, the imaginary part e" of its dielectric 

function e(A) = e' + i e" is expected to be larger compared to that 
of its bulk counterpart, as a consequence of the so-called electron 
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mean free path limitation effect (e.g. see Boliren & Huffman 1983). 
Tliis is easier to understand if we decompose e into two components 
66 and ef, contributed by bound charges ("interband transitions") 
and free electrons, respectively. The free electron component is well 

described by the Drude theory e = 1 — ^•2_^^- ^ with an imaginary 

2 

part e" = ^^.J^^i^^ where the plasma frequency cOp is related to the 

free electron density Ue- cOp = neC^/m*, the damping constant 7 is 
related to the average time r between collisions: 7=l/r. In bulk ma- 
terials, 7 is mainly determined by the scattering of the electrons with 
phonons (lattice vibrations), and to a lesser degree, with electrons, 
lattice defects, or impurities. However, for particles in the nanome- 
ter size domain, 7 is increased because of the additional collisions of 
the conducting electrons with the grain boundary: 7 = 7buik + 
where 7bulk is the bulk metal damping constant, vp is the electron 
velocity at the Fermi surface, /3 is a dimensionless constant of order 
unity which depends on the character of the scattering at the bound- 
ary: /3=1 for classic isotropic scattering, (3=4:/ 3 for classic diffusive 
scattering, /3=1.16 or /3=1.33 for scattering based on the quantum 
particle-in-a-box model (see Coronado &z Schatz 2003 and references 
therein). Since uj"^ » 7^ in metals near the plasma frequency, e" can 

2 

be written as e" = e'-^^^^ + ^^^g . Clearly, for a metallic grain e" in- 
creases as the grain becomes smaller. 

In contrast, our knowledge regarding the size dependence of the di- 
electric function of dielectric materials is controversial. For example, 
some theoretical and experimental studies have concluded that the 
dielectric function of Si nanoparticles is significantly reduced relative 
to the bulk value (Koshida et al. 1993; Wang & Zunger 1994; Tsu, 
Babic, & loriatti 1997; Amans et al. 2003), which is generally at- 
tributed to the quantum confinement effect. But we find that the 
absorption and reflectivity measurements for porous silicon appear 
to be consistent with bulk Si together with voids and Si02 (see Li 8z 
Drainc 2002a for details). 

— It has been reported that the specific heat of small metal particles is 
strongly enhanced over the bulk value (see Halperin 1986, Meyer et al. 
2003 and references therein). For example, a progressive decreasing 
of Debye temperature G with the decrease of grain size has been 
observed for palladium (Pd): 273, 226, 193, 175 K for bulk Pd and 
Pd particles of radius a = 42, 33, 15 A, respectively (Chen et al. 1995 
and references therein). This has been attributed to quantum effects 
on the vibrational spectrum in small particles: as a grain becomes 
smaller, a larger fraction of atoms occupy surface sites which are 
weakly bound to the grain; therefore, small grains are expected to 
have a larger low-frequency mode density due to the weaker bonds of 
the surface atoms. 

However, although an enhancement of the vibrational specific heat 
is also expected for dielectric grains, the degree of enhancement is 
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unclear due to the differences between the binding and structural 
properties of dielectrics with metals. For example, based on a lattice 
dynamical calculation, Hu et al. (2001) have shown that the differ- 
ence in vibrational specific heat between Si nanocrystals and the bulk 
is just about a few per cent. The specific heat of nanocrystalline di- 
amond measured by Moelle et al. (1998) shows a close agreement 
with that of hulk diamond. It is worth noting that the heat capac- 
ity of PAHs can be well described by the two 2-dimcnsional Dcbye 
models of hulk graphite (together with contributions from the C-H 
vibrational modes) (see Fig. 2 of Draine & Li 2001). Clearly, lab- 
oratory studies of the optical and thermal properties of interstellar 
nanoparticlc analogues would be very valuable. 
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